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Abstract— We analyze the encounters of the neutron star (pulsar) Geminga with open 
star clusters in the OB association OriOBla through the integration of epicyclic orbits into 
the past by taking into account the errors in the data. The open cluster ASCC21 is shown to 
be the most probable birthplace of either a single progenitor star for the Geminga pulsar or a 
binary progenitor system that subsequently broke up. Monte Carlo simulations of Geminga- 
ASCC21 encounters with the pulsar radial velocity Vr = —100 ± 50 km s~^ have shown that 
close encounters could occur between them within < 10 pc at about t = —0.52 Myr. In 
addition, the trajectory of the neutron star Geminga passes at a distance of ~25 pc from 
the center of the compact OB association A Ori at about t = —0.39 Myr, which is close to 
the age of the pulsar estimated from its timing. 

INTRODUCTION 

Young massive O and B stars can run away from the stellar group (an open cluster or an 
OB association) where they were born. Such stars are called runaways; they are identified 
by high (> 30 km s~^) space velocities. Two main runaway scenarios are known (see the 
review by Hoogerwerf et al. 2001). The first scenario is related to the evolution of a binary 
system and the explosion of one of the components as a supernova. In this case, one of the 
components can receive a high kick velocity due to an asymmetric explosion and a significant 
mass loss by the binary. A massive runaway star and its former companion, a neutron star 
or a black hole, are the breakup products of such a binary. The second scenario is related to 
the dynamical evolution of an initial stellar group. Mutual encounters of close binary and 
multiple systems are the most efficient ejection mechanism. 

The best-known runaway stars, such as COph, ^Per, /xCol, AE Aur, CPup, and others, 
were discovered by Blaauw (1961; 1993). Probable parent associations were determined for 
some of them. Hoogerwerf et al. (2001) considered the problem of determining the most 
probable birthplaces of pulsars using data from the Hipparcos (1997) Catalogue. Based on a 
more complete list of runaway OB stars, Schilbach and Roser (2008) showed that the parent 
open cluster or OB association to which they belonged could be pinpointed with a high 
probability for most (~ 90%) of the O stars in the solar neighborhood (within < 1 kpc), 
including those that were previously believed to be field stars. 

The problem of determining probable birthplaces of nearby young neutron stars associ- 
ated with the Gould Belt (Popov et al. 2003; Motch et al. 2006) is of great interest. One of 
these stars is Geminga (PSR J0633+1746), to which this paper is devoted. 



In contrast to the overwhelming majority of pulsars, Geminga was discovered in 1973 
as a gamma-ray source from the SAS-2 sateUite. The impulsive behavior of its signals was 
first detected in the X-ray (Halpern and Holt 1992) and gamma-ray (Burch et al. 1992; 
Bignami and Caraveo 1992) bands and, only several years later, in the radio band (Kuzmin 
and Losovsky 1997; Malofeev and Malov 1997). The age of the Geminga pulsar, 0.342 Myr, 
was first determined from its timing (P/2P , where P and P are the pulsar period and its first 
derivative, respectively) by Bignami and Caraveo (1992). Bignami et al. (1993) identified 
this pulsar in the optical band with a 25.5-magnitude object and were the first to determine 
the components of its proper motion. The distance to the Geminga pulsar, r = 150 — 400 
pc, was first estimated by Halpern and Ruderman (1993). 

The present-day quality of the spectra for pulsars is too low for their radial velocities to 
be determined. At the same time, an arc-like feature produced by the interaction of pulsar 
emission with the interstellar medium is observed near the neutron star Geminga (Caraveo 
et al. 2003; Shibanov et al. 2006). Modeling the shape of this arc-like feature (Caraveo 
et al. 2003) yielded a constraint on the magnitude of the radial velocity for the Geminga 
pulsar {\Vr\ < 100 km s~^). Given this constraint, the radial velocity is a free parameter in 
our studies. 

One of the first hypotheses that the birthplace of the Geminga pulsar is associated with 
the Local Bubble was proposed by Gehrels and Chen (1993). According to this hypothesis, a 
positive radial velocity of the pulsar should be adopted. In the opinion of Frisch (1993), this 
pulsar originated in the Orion OB association, suggesting that its radial velocity is negative. 
The modeling by Smith et al. (1994) showed that the compact OB association AOri (Cr 69) 
at K ~ —100 km s~^ could be a suitable birthplace of the pulsar. 

Using Hubble Space Telescope (HST) observations, Caraveo ct al. (1996) were the first 
to determine the trigonometric parallax of the pulsar, TTtr = 6.36 ± 1.74 mas, which corre- 
sponds to its distance r = 1571^4 PC- They also determined the proper motion components, 
Ha cos S — 138±4 mas yr~^ and hs = 97±4 mas yr~^. The search for a probable birthplace of 
the neutron star Geminga carried out by PeUizza et al. (2005) based on these data led them 
to conclude that this birthplace could be in a fairly close solar neighborhood (r = 90 — 240 
pc) and is related either to the relatively old (^50 Myr) CasTau association or to the edge 
of the younger (aill Myr) OB association Ori OBla. 

The currently available HST observations of Geminga (Faherty et al. 2007) showed 
that the new parallax, ntr — 4.0 ± 1.3 mas (r = 250t.ll° pc), differs significantly from the 
previously published one (Caraveo et al. 1996), which should inevitably affect the results of 
searching for its birthplace. 

The goal of this paper is to determine the most probable birthplace of the Geminga 
pulsar by simulating its motion in space into the past using the most resent data on open 
star clusters (Kharchenko et al. 2007) and individual stars (van Leeuwen 2007). The objec- 
tives set here are accomplished using the method of statistical Monte Carlo simulations-we 
determine the epicyclic orbits of objects by taking into account the corridor of errors in the 
observational data. 

DATA 



The input data on the objects under consideration, such as their equatorial coordinates, 
proper motion components, radial velocities, and parallaxes, are given in the table. The data 



Table 1: Data on stars and open star clusters Object. 



Object 


Q;(J2000.0); 


fla COS 5, 






Distance, 




^(J2000.0) 


fj,s. mas yr^"^ 


mas 


km s^^ 


pc 


ASCC16 


5/i24m3gs 
1°48'00" 


0.75 ±0.22 
-0.18 ±0.29 




20.80 ±4.60 


460 
(2) 


ASCC18 


5''26'"09.6'^ 
0°49'12" 


0.89 ±0.28 
-0.02 ±0.25 




24± 10 


500 

(2) 


ASCC20 


1°37'48" 


-0.09 ±0.21 
0.51 ±0.19 




22.97 ±4.80 


450 

(2) 


ASCC21 


5''28™58.8" 
3°39'00" 


0.52 ±0.25 
-0.62 ±0.28 




19.77 ± 1.12 


500 
(2) 


A Ori 


5'*35'"09.6* 


0.40 ± 0.47 




31.38 ± 1.42 


438 


(Cr 69) 


9°42'00" 


-1.93 ±0.34 






(2) 


HIP 22061 


4M4'"42M571 


-44.02 ± 0.52 


2.99 ±0.63 


9.0 ±4.4 






0°34'05".418 


-30.04 ±0.50 


(1) 


(3) 




HIP 29678 


6''15'^08^4567 


25.23 ±0.50 


2.72 ± 0.43 


35.8 ± 1.0 






13°51'03".859 


11.44 ±0.30 


(1) 


(4) 




Geminga 


6'^33"'54M530 
17°46'12".909 


142.2 ± 1.2 
107.4 ± 1.2 


4.0 ± 1.3 

(5) 







Note. (1) van Leeuwen (2008), (2) CRVOCA (Kharchenko et ak 2007), (3) PCRV (Gontcharov 
2006), (4) CRVAD-2 (Kharchenko et ak 2007), (5) Faherty et al. (2007). 



for all of the stars used were taken from a revised version of the Hipparcos Catalogue (van 
Leeuwen 2007). For the Geminga pulsar, wc use new data from Faherty et al. (2007), which 
were obtained from 18-month-long HST observations of the star. The CRVOCA Catalog 
(Kharchenko et al. 2007) served as the source of data for open clusters. The distance 
estimates in this catalog are based on a photometric method with an error of ^ 20%. We 
calculated the parallaxes for each open cluster based on the distances given in the last column 
of the table. 

Four groups, Ori OBla, b, c, and d (Blaauwl964), with ages of 11.4±1.9, 1.7±1.1, 4.6±2, 
and < 1 Myr, respectively, are identified in the Orion association (Brown et al. 1994). 
According to Kharchenko et al. (2005), the Ori OBla association includes four open clusters, 
ASGG16, ASGC18, ASGG20, and ASGG21, with ages of 8.5, 13.2, 22.4, and 12.9 Myr, 
respectively. The question of whether each of them is an open cluster or an association is 
still under discussion. In this respect, ASCC16, which is also designated as the group 25 Ori 
(Briceno et al. 2005; McGehee 2006), is the best-studied open cluster. Here, a fairly dense 
concentration of young T Tauri stars was detected. Their analysis led McGehee (2006) to 
conclude that this is a gravitationally unbound group rather than an open cluster. 

The OB-T association/cluster AOri is located near the Ori OBla association. A compact 
HII region (S264) is associated with the star AOri. In a wider neighborhood of the center 
with a radius of ~ 5° (40 pc), the association is surrounded by a ring of molecular clouds and 
gas (Maddalena and Morris 1987). This star- forming region was studied in a series of papers 
by Dolan and Mathieu (1999, 2001, 2002). Based on Stromgren photometry, they estimated 



the distance to the open cluster AOri, 450 ± 50 pc. Ten B stars forming the cluster proper 
were shown to concentrate in a sky field with a radius of ~ 1° (8 pc) around the central star 
of the cluster AOri (Sp: 08III) with an age of Ri5.5 Myr. In a wider neighborhood of the 
center with a radius of ~ 5° (40 pc), they detected several concentrations of young T Tauri 
stars. On the whole, Dolan and Mathhieu concluded that: (a) star formation in the cluster 
AOri ceased ~1 Myr ago after the supernova explosion; (b) the group of young massive stars 
is currently not gravitationally bound; and (c) star formation currently goes on in the clouds 
B30 and B35 located on the periphery of the association. 



METHODS 

Orbit Construction 

In this paper, we use a rectangular Galactic coordinate system with the axes directed away 
from the observer toward the Galactic center (/ = 0°, 6 = 0°, the X axis), along the Galactic 
rotation (/ = 90°, 6 = 0°, the Y axis), and toward the North Galactic Pole (h = 90°, the Z 
axis). The corresponding space velocity components of the objects U, V, W are also directed 
along the X, Y, Z axes (Kulikovskii 1985). 

The epicychc approximation (Lindblad 1927, 1959) allows the orbits of stars to be con- 
structed in a coordinate system rotating around the Galactic center. We assume that the 
origin of the coordinate system coincides with the local standard of rest and that the stars 
move along epicycles in the direction opposite to the Galactic rotation. We use the method 
in the form given by Fuchs et al. (2006): 

X(t) = X(0) + ^ sm(Kt) + ^(1 - cos{Kt)), (1) 

Y{t) = y(0) + 2a(x{0) + - ^V{0) sm{Kt) + ^C/(0)(1 - cos{Kt)), 

Z{t) = y^ sin(i/t) + Z(0) cos(i/t), 

where t is the time in Myr (pc/Myr=0.978 km s~^), which is measured into the past; A 
and B are the Oort constants; k — ^/—4:QoB is the epicyclic frequency; Qq is the angular 
velocity of Galactic rotation of the local standard of rest, Qo = A — B; and u = y^ArcGpo is 
the frequency of the vertical oscillations, where G is the gravitational constant and po is the 
star density in the solar neighborhood. The space velocities of the objects are calculated for 
any necessary time using the formulas 

U{t) = U{0) cos{Kt) + :^V{0) sm{Kt) , (2) 

2x3 

Vit) = 2^(X(0) + ^) - ^ViO) cos(«:t) + ^^(0) sin(«:t)), 
W{t) = "^^(0) cos(i/t) - Z(0)i/sin(i/i). 

The parameters X{0),Y{0), Z{0) and U{0),V{0),W{0) in Eqs. (1) and (2) denote, respec- 
tively, the current positions and velocities of the objects. The velocities U, V, W are given 
relative to the local standard of rest with {U, V, W)lsr = (10.00, 5.25, 7.17) km s^-^ (Dehnen 
and Binney 1998). Following Fuchs et al. (2006), we adopted po = O.IM© pc^, which gives 



1/ = 74 km kpc~^. We also used the Oort constants A = 13.7 ± 0.6 km s^^ kpc^^ and 
B — —12.9 ± 0.4 km kpc~^ that we found previously (Bobylev 2004) by analyzing the 
independent determinations of these parameters by various authors; k — 37 km s~^ kpc~^ 
corresponds to these values. 



Statistical Simulations 

In accordance with the method of statistical Monte Carlo simulations, we compute 3 million 
orbits for each object by taking into account the random errors in the input data. For each 
pair of orbits belonging to two different objects, we calculate the encounter parameter equal 



to the minimum separation between the objects Ar = AX'^(t) + AV^it) + AZ'^(t). 

The parameters of the objects are distributed normally with a dispersion a. As was 
shown by Hoogerwerf et al. (2001), the random errors in the equatorial coordinates of stars 
and open clusters do not affect noticeably the simulation results. As a result, the errors 
are added only to the proper-motion components, parallax, and radial velocity of an object. 
The parameters are listed in the table. 

The distribution of the encounter parameter in a certain neighborhood can be represented 
as a histogram (see below). We calculated the expected distribution F^d of the minimum 
separation using a formula from Hoogerwerf et al. (2001), 



for appropriate mean /j, and dispersion a. When /i may be considered close to zero, Hooger- 
werf et al. (2001) obtained the expression 



2 



A 



4a2 



(4) 



For such open clusters as ASCC21 and AOri, the tidal radius is 13.7 ± 1.8 and 7.7 ± 1.4 pc, 
respectively (Piskunov et al. 2008). Therefore, in our simulations, we choose the character- 
istic radius of the neighborhood in which close encounters of stars with clusters occur to be 
10 pc. 

Analysis of the orientation of the shock front ahead of the pulsar yields a constraint on 
its radial velocity. A fairly symmetric shock was detected near the neutron star Geminga 
(Caraveo et al. 2003; Shibanov et al. 2006). Modeling its shape suggests that the radial 
velocity of Geminga is relatively low, l^^l < 72 km s"-*^; the angle /3 between the total space 
(V) velocity vector and the tangential {V±) velocity is < 30° (Garaveo et al. 2003; PeUizza 
et al. 2005). For the new data (table), f3 = 30°. 7, where |cos/3| = V±/V, corresponds to 
the radial velocity = 100 km s~^. Therefore, for the Geminga pulsar, we take its model 
radial velocity Vr = —100 ± 50 km s~^. 



RESULTS 

Figure 1 shows the positions and trajectories of the open clusters in the Ori OBla association, 
the cluster AOri, and the neutron star Geminga for a radial velocity of —100 ± 50 km s~^ 
over 1 Myr into the past (with marks at 0.5-Myr intervals). To construct the trajectories 
in Galactic coordinates, we first calculated the rectangular coordinates X{t), Y{t), and Z{t) 



from Eqs. (1) and then found the corresponding spherical Galactic coordinates l{t) and b{t) 
from them. As can be seen from the figure, the trajectory of the Geminga pulsar passes near 
two open clusters, ASCC21 and AOri. The two dotted lines in Fig. 1 indicate the pulsar 
trajectories calculated with extreme (±3cr) values of the proper-motion components. 

Our simulations of the Geminga encounter with the cluster ASCC21 yielded the following 
results. Out of the 3 million orbits obtained by introducing random errors with a given 
dispersion for both pulsar and cluster, 10 5441 encounters occur at separations < 10 pc 
(3.5%). In 3137 of the 10 5441 cases (3.0%), the encounters at about t = -0.52 Myr occur 
in the neighborhood of ASCC21 with a radius < 10 pc and with the center lying on the 
cluster trajectory obtained from the data without introducing any errors. Figure 2a shows 
the expected distribution F^d of the minimum separation calculated from Eq. (3) for the 
adopted mean /i — 7.75 pc and dispersion cr = 1.0 pc. The domain of admissible values of V^, 
TT, /la cos 5 {iJ>a), and /is at which 3137 encounters occur are shown in Fig. 3 for the Geminga 
pulsar and Fig. 4 for the cluster ASCC21. The horizontal and vertical fines in Figs. 3 and 4 
pass through the nominal values of the parameters used. 

Figure 3 shows the set of lines representing the region that satisfies the condition < 30° 
for the Vr—TT relation. The fines were obtained in accordance with the relation K = Vl. tan (3. 
The solid and dotted lines correspond to the angles P — ±30° and /3 = ±15°, respectively. 
As we see from the figure, most of the model points fall into the interval < 30°, thereby 
satisfying the constraint on the pulsar radial velocity determined from an analysis of the 
shock shape (Caraveo et al. 2003; Pellizza et al. 2005). 

Our simulations of the encounter of the neutron star Geminga with the cluster AOri re- 
vealed no mutual encounters at separations <10 pc. However, there are less close encounters. 
Thus, for example, out of the 300 000 model orbits, 25 745 encounters occur at separations 
Ar < 25 pc (8.5%); in 2686 of the 25 745 cases, the encounters occur within <30 pc (10%) 
of the cluster AOri at about t = —0.39 Myr. 

Our simulations of the Geminga encounter with the cluster ASCC16 yielded the following 
results: out of the 3 million orbits, 27 848 encounters occur at separations A^ < 10 pc (0.9%); 
in 1475 of the 27 848 cases, the encounters occur within <50 pc of the cluster ASCC16 at 
about t = —0.63 Myr. Thus, in this case, we have a factor of 4 fewer encounters than in the 
case of ASCC21. 

To determine a possible candidate for the binary component with the neutron star 
Geminga in the past, we selected known rapidly flying stars whose trajectories passed near 
the Ori OBI association on the celestial sphere in a time interval of 2 Myr (Fig. 2 from 
Hoogerwerf et al. 2001): HIP 22061, HIP 27204, and HIP 29678 (numbers 4, 6, and 7 in the 
list of Hoogerwerf et al., 2001). 

For HIP 22061, there are encounters with the Geminga pulsar at i 0.6 Myr but within 
Ar < 20 pc, i.e., the encounters are not close. For HIP 27204, there are no encounters with 
the Geminga pulsar within A,. < 20 pc. For HIP 29678, there are close encounters with the 
Geminga pulsar. We found that out of the 3 million orbits, 124 544 encounters occurred at 
separations A^ < 10 pc (4%) at t = —0.14 Myr, which is half the pulsar age calculated from 
its timing. 

According to Hoogerwerf et al. (2001), the stars HIP 22061 (B2.5V) and HIP 29678 

(BIV) have masses of 8.6M0 and 11.5AfQ, respectively. They are probably the components 
of a dissociated multiple (HIP 29678 is a visual binary) system. We found that out of the 
300 000 orbits, 18 165 encounters occur between these two stars at separations A^ < 10 
pc (6%); in addition, 5310 of the 18165 encounters occur within A^ < 50 pc of the cluster 



AOri at about t = —1.12 Myr. As can be seen from the table, the locahzation error for 
each of the objects under consideration is 20%. This served as a basis for choosing the 
size of the neighborhood (50 pc). Therefore, our results are consistent with the assumption 
that HIP 22061 and HIP 29678 ran away from the neighborhood of the cluster AOri about 
1.1 Myr ago. Figure 5a shows the expected distribution F3D of the minimum separation 
Ar calculated from Eq. (4) for the dispersion a = 1.2 pc. This is an example of very close 
encounters; thus, for example, the separation between the two stars was only <0.12 pc in 
two model cases. 

The conclusion that HIP 22061 and HIP 29678 are the remnants of a dissociated binary or 
multiple system is in agreement with the analysis of the motion of these stars performed by 
Hoogerwerf et al. (2001) using Hipparcos data. However, Hoogerwerf et al. (2001) questioned 
the involvement of the star AOri in the event that led to the breakup of the multiple system 
HIP 22061-HIP 29678. This conclusion was drawn from the fact that the point on the 
celestial sphere with coordinates {l,b) = (196°. 5, —12°. 0) was found as a probable current 
position of the parent cluster for HIP 22061 and HIP 29678, which differs by one degree from 
the current position of the center of the cluster AOri. The velocity difference effect can be 
seen in Fig. 1. More specifically, if the cluster AOri were located one degree to the left, then 
all three trajectories (AOri, HIP 22061, and HIP 29678) would intersect at the same time, 
t « -1 Myr. 

The kinematic paradox in the motion of the star AOri was pointed out by Maddalena 
and Morris (1987). It lies in the fact that the position of the star AOri in the past does not 
coincide with the center of the ring of molecular hydrogen clouds expanding with a velocity 
of 14.3±2.5 km s~^. We are talking about a velocity difference of 7 — 10 km s"^. According to 
Maddalena and Morris (1987), the coordinates of the ring center are {a, S)ig5o — (5'*29"*.8± 
1".6,9°54' ± 24') and, hence, {l,b) = (194°.7 ± 0°.4, -12°.5 ± 0°.4) (in our view, these 
coordinates reflect the position of the supernova at the time of its explosion). The position 
of the ring center remarkably coincides with the point of intersection of the trajectories for 
HIP 22061 and HIP 29678 that we found, namely, {l,b) = (194°.6 ± 0°.2, -12°.2 ± 0°.2). 
The errors here were estimated from the relation = arctan((7/r), where r = 438 pc, and 
from the dispersion o" = 1.2 pc found above (Fig. 5a). 

Let us assume that the HIP 22061-HIP 29678 pair is a relic of the system (the cluster 
AOri) before the time that preceded the dynamical "event" as a result of which the star AOri 
received an additional kick (—10 km s~^). In this case, we obtain an independent estimate 
at the time of the event — it could occur no later than 1.12 Myr ago. Moreover, it follows 
from the above assumption that the point of intersection of the HIP 22061-HIP 29678 pair 
indicates the position of the center of the cluster AOri at a time of —1.12 Myr. This suggests 
that the supernova explosion that gave rise to the ring of molecular clouds around the cluster 
AOri occurred no later than —1.12 Myr. 

We can see that the pair of runaway starsHIP22061 and HIP 29678 is of great interest in 
studying the evolution of the cluster AOri. We consider the close encounters of HIP 29678 
with the Geminga pulsar found above as a purely geometric effect. If both these objects 
(HIP 29678 and Geminga) resulted from the evolution of a binary progenitor system, then 
its trajectory before its breakup could be arbitrary and the encounter of the HIP 22061-HIP 
29678 trajectories would not be as shown in Fig. 1. 



1 DISCUSSION 



Based on an isochronous age estimate for the cluster ASCC21, 12.9 Myr (Kharchenko et al. 
2005), we can estimate the minimum mass of a post-main sequence star. For example, based 
on Fig. 11.5 from SakhibuUin (2003) for logt = 7.11, we obtain an estimate of the mass for 
such a star, ^ 15Mq. On the other hand, there are no stars bluer than a B1.5V-type star 
(HIP 25861) in the cluster ASCC21. Therefore, the more massive BOV-type stars, whose 
masses are known to be ^ 15M0, can be assumed to have already left the main sequence. 
A more detailed comparison with isochrones, for example, from Schaller et al (1992), also 
leads to a similar conclusion. 

We may conclude that in the case of ASCC21, the mass of the probable progenitor star 
for the Geminga pulsar satisfies well the condition for the formation of neutron stars with 
(10 - 12)M0 < M < (30 - 4O)M0 (Zasov and Postnov 2006). 

The cluster AOri is considerably younger. Dolan and Mathieu (2002) estimated the mass 
and age of the most massive cluster star, AOri, to be 26. SM© and a; 5.5 Myr, respectively. 
Thus, the higher-mass stars have already evolved. Therefore, the cluster ASCC21 appears 
more preferable in the scenario with a single progenitor star for the Geminga pulsar. 

Our simulations showed that there is a probability (~ 3%) of the simple scenario in 
which the pulsar was formed in ASCC21 aX t — —0.52 Myr and has since moved in the same 
direction with the initially received velocity. 

As was pointed out by Yakovlev et al. (1999), the characteristic (i.e., dynamical, de- 
termined from timing) age of Geminga was estimated to within a factor of ~3. Glitches in 
period, which are also observed in the Geminga pulsar (Jackson et al. 2002), are among the 
factors that affect the dynamical age estimates for pulsars. On the other hand, the neutron 
star Geminga exhibits thermal radiation, which allows its age boundaries to be estimated 
by comparison with neutron star cooling models. Thus, for example, we can see from the 
data in Table 3 and the cooling curves in Fig. 20 from Yakovlev et al. (1999) that the upper 
limit for the Geminga age in the blackbody model reaches ~ 1 Myr. Similar estimates of 
the age boundaries for the Geminga pulsar, 0.1 — 1 Myr, can be found in Page et al. (2004). 

Thus, the scenario for the formation of Geminga in the cluster ASCC21 at t = —0.52 
Myr is consistent with the time intervals considered. For this scenario to be realized, first, 
the Geminga parallax must fall within the range 2-3 mas (Fig. 3) and, second, there must 
be traces of the supernova remnant around ASGG21, which is not yet observed. On the 
other hand, such a remnant could disperse completely in a 0.5 Myr. 

The assumption that the progenitor star from which the Geminga pulsar subsequently 
originated was formed in ASCC21 appears more realistic. This viewpoint is in agreement 
with the scenario by Smith et al. (1994), who believe the Geminga pulsar to be either the 
remnant of a runaway OB star ejected from the Orion OBla association and exploded as 
a supernova outside this group or the result of a supernova explosion in the cluster AOri. 
Cunha and Smith (1996) provide arguments that the ring of molecular hydrogen clouds 
around the cluster AOri has an age of 300-370 thousand years, in good agreement with the 
dynamical age estimate for the neutron star Geminga. 



CONCLUSIONS 



Analysis of the possible encounters of the neutron star Geminga with open clusters in the Ori 
OBla association in the past led us to conclude that the open cluster ASCC21. a member of 
the Ori OBla association, is the most likely candidate for the birthplace of either a binary 
progenitor system that subsequently broke up or a single progenitor star. Monte Carlo 
simulations of the encounters of this pair with Gemingas radial velocity K = —100 ± 50 
km s~^ showed that close encounters of the pulsar with this open cluster could occur within 
< 10 pc at about t = —0.52 Myr. In addition, the trajectory of the Geminga pulsar passes 
at a distance of about 25 pc from the center of the compact OB association AOri at about 
t = -0.39 Myr. 

Previously, the Geminga pulsar was assumed to be related to a wide neighborhood of the 
Ori OBla association, while we found the specific open cluster (ASCC21) from which the 
massive progenitor star or the progenitor binary could be ejected. 

In addition, the kinematics of the pair of runaway stars HIP 22061-HIP 29678 under the 
assumption that they were members of the cluster AOri in the past allowed us to obtain an 
independent estimate for the time of the dynamical "event" as a result of which the star 
AOri received an additional kick — it could occur no later than 1.12 Myr ago. 
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Fig. 1. Positions of the open clusters in the Ori OBI association and the Geminga pulsar 
(the trajectory was computed for a radial velocity of —100 km s~^). 
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Fig. 2. (a) Expected distribution of the minimum separation A,. < 10 pc for Geminga 
encounters with the cluster ASGG21; (b) histogram of encounter times (the time is measured 
into the past). 
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Fig. 3. Domains of admissible values for Geminga at which 3137 encounters occur between 
the pulsar and the cluster ASCC21 to distances < 10 pc. 




Fig. 4. Domains of admissible values for the cluster ASCC21 at which 3137 encounters with 
the pulsar occur at separations < 10 pc. 
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Fig. 5. (a) Expected distribution of the minimum separation A,. < 10 pc for the encounters of 
HIP 22061 and HIP 29678 within < 50 pc of the cluster AOri; (b) histogram of encounter 
times (the time is measured into the past). 



